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We present a nonparametric approach to reconstruct the interaction between dark energy and 
dark matter directly from SNIa Union 2.1 data using Gaussian processes, which is a fully Bayesian 
approach for smoothing data. In this method, once the equation of state {w) of dark energy is 
specified, the interaction can be reconstructed as a function of redshift. For the decaying vacuum 
energy case with w = —1, the reconstructed interaction is consistent with the standard ACDM 
model, namely, there is no evidence for the interaction. This also holds for the constant w cases from 
—0.9 to —1.1 and for the Chevallier-Polarski-Linder (GPL) parametrization case. If the equation of 
state deviates obviously from —1, the reconstructed interaction exists at 95% confidence level. This 
shows the degeneracy between the interaction and the equation of state of dark energy when they 
get constraints from the observational data. 

PACS numbers: 95.36.+X, 98.80.-k, 98.80.Es 


I. INTRODUCTION 

It has been more than fifteen years since the universe 
was found in accelerating expansion [U, 0. However, it is 
fair to say that its origin is still not yet clear. A possi¬ 
ble explanation of this cosmic acceleration is provided by 
the introduction of a fluid with negative pressure called 
dark energy (DE). The simplest dark energy candidate 
is the cosmological constant A with the equation of state 
w = —1. The tiny cosmological constant together with 
the cold dark matter (CDM) (called the ACDM model) 
turned out to be the standard model which fits the cur¬ 
rent observational data sets consistently. In spite of this 
success, however, it is faced with the fine-tuning prob¬ 
lem Q and the coincidence problem. The former arises 
from the fact that the present-time observed value for 
the vacuum energy density is more than 120 orders of 
magnitude smaller than the naive estimate from quan¬ 
tum field theory. The later is the question why we live in 
such a special moment that the densities of dark energy 
and dark matter are of the same order. 

Many attempts have been made to tackle those issues, 
including introducing “dynamical” dark energy or mod¬ 
ifying general relativity at the cosmic scales. In particu¬ 
lar, to alleviate the coincidence problem, an interaction 
between dark energy and dark matter has been consid¬ 
ered extensively in the literature, for example, in 1M3 
and references therein. However, usually the interaction 
form has to be assumed. The assumed form will lead 
to a bias when observational data are used to give the 
constraint on the interaction. In the authors inves¬ 
tigated the possible interaction in a way independent of 
specific interacting forms by dividing the whole range of 
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redshift into a few bins and setting the interacting term to 
be a constant in each redshift bin. It was found that the 
interaction is likely to cross the noninteracting line and 
has an oscillation behavior. Recently, Salvatelli et al. [Tsjl 
showed that the null interaction is excluded at 99% con¬ 
fidence level (C.L.) when they added the redshift-space 
distortions (RSD) data to the Planck data for the decay¬ 
ing vacuum energy model (a class of interaction of dark 
sectors). They parametrized the interaction term to be 
of the form qHp and also subdivided the redshift into 
four bins with q(z) = qi(i = 1,... ,A). More recently, the 
authors of Ref. [I^ have reconstructed the temporal evo¬ 
lution of the coupling strength between dark matter and 
vacuum energy, a (a) in a nonparametric Bayesian ap¬ 
proach using the combined observational data sets from 
the cosmic microwave background radiation, supernovae 
and large scale structure. It was found that an evolving 
interaction can remove some of the tensions between dif¬ 
ferent types of data sets, and is favored at ~ 95% G.L. 
if the baryon acoustic oscillations measurements of the 
BOSS Lyman-a forest sample are included. Thus, it is 
quite interesting to see whether there exists some signa¬ 
ture of the interaction from the observational data in a 
model-independent way. 

In this paper we present a nonparametric approach 
to reconstruct the interaction term between dark energy 
and dark matter directly from the observational data us¬ 
ing Gaussian processes (GP). GP is a model independent 
method to smooth the data. We set the nonparameter- 
ized interaction term Q{z) as a function of redshift and 
reconstruct it from SNIa Union 2.1 data sets. We will 
consider three cases, the decaying vacuum energy case 
with w = —1, the rcGDM model and the Ghevallier- 
Polarski-Linder (GPL) parametrization of dark energy, 
respectively, and pay special attention to the first case as 
discussed in [l5l| . 

This paper is organized as follows. In Sec. |TT1 "we give 
the interacting models of dark energy and dark matter 
in a flat universe. In Sec. uni we briefly introduce the 
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Gaussian processes method and create a mock data set 
to demonstrate the reliability of the GP reconstruction 
method. Then we apply it to the real data Union 2.1 in 
the decaying vacuum energy case, followed by different 
values of w for comparison and the GPL case. We give 
some discussions and make conclusions in Sec. m 


II. THE MODEL 

In a flat universe with an interaction between dark en¬ 
ergy and dark matter, the Friedmann equation describing 
the evolution of the universe is given by 

r, STTG , . , . 

H = —^{Pm+PDE), ( 1 ) 


The luminosity distances at redshift 2 : can be expressed 
as 
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Writing D{z) = {Ho/c){l -I- z) ^dhiz) as the normalized 
comoving distance, we can arrive at 
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where pm denotes the energy density of dark matter and 
Pde the energy density of dark energy. However, the 
conservation equations are changed to be 

Pm + ^Hpm = —Q, (2) 


From this, we see that using the observed distance- 
redshift relationship D{z), one can reconstruct the in¬ 
teraction, once the equation of state w of dark energy is 
given. 


Pde+ '^H{1 + w)pde = Q, (3) 


III. RECONSTRUCTION METHOD 


where w is the equation of state of DE, H is the expan¬ 
sion rate of the universe, and Q describes the interaction 
between dark matter and dark energy. When (5 = 0 and 
w = — 1, it recovers the standard ACDM model. Unlike 
most of the ways to parametrize the interaction term Q 
using such a form qHp [l,[iH, here we use Q{z) itself and 
want to reconstruct it directly from data using a model- 
independent method. As assumed in [l^, the 4-vector Q'" 
is proportional to the 4-velocity of dark matter. Hence 
the perturbation of the interaction is not considered in 
this paper. 

Combining the Friedmann equation and the conserva¬ 
tion equations {SttG = 1 throughout the paper), we can 
obtain 


-wQ = 2 ( HH'^ + H^H" - —H^H' ) (1 -F zf 
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where the prime denotes the derivative with respect to 
redshift z. Note that here we also assume w is time de¬ 
pendent. For convenience, we use a dimensionless q to 
characterize the interaction, i.e., Q = qH^. Note that 
here q is not the usual deceleration. In this case, we have 


—wqHo^ = 2 


-h H^H" 


- 2{5+ 3w)H^H' 
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In order to reconstruct the interaction using current 
data sets, we should find a mo del-independent method 
to reconstruct D{z) and its derivatives. While there 
are several methods such as principle c omp onent anal¬ 
ysis (l9l - [^ . Gaussian smoothing and Gaus¬ 

sian processes l2ii3> in this paper we will reconstruct 
D{z) and its derivatives more precisely by using the GP 
method. 


A. Gaussian processes 

The Gaussian processes allows one to reconstruct a 
function from data without assuming a parametrization 
for it. We use Gaussian processes in Python (GaPP) [l^ 
to derive our GP results. This GP code has been applied 
in many papers (27l - [^ . The distribution over functions 
provided by GP is suitable to describe the observed data. 
At each point z, the reconstructed function /(z) is also 
a Gaussian distribution with a mean value and Gaus¬ 
sian error. The functions at different points z and z are 
related by a covariance function fc(z, z), which only de¬ 
pends on a set of hyperparameters i and af. Here i gives 
a measure of the coherence length of the correlation in 
x-direction and cr/ denotes the overall amplitude of the 
correlation in the y-direction. Both of them will be op¬ 
timized by GP with the observed data set. In contrast 
to actual parameters, GP does not specify the form of 
the reconstructed function. Instead it characterizes the 
typical changes of the function. 

The different choices of the covariance function may 
affect the reconstruction to some extent. The covariance 
function is usually adopted as the squared exponential 
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form 

fc(z, 5) = exp ( - . (8) 

But it is not always the suitable choice. Here we take the 
Matern {v = 9/2) covariance function 


k{z, z) = exp I — 
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according to the analysis made in (^ . where they con¬ 
sider various assumed models and many realizations of 
mock data sets for a test and conclude that the Matern 
{v = 9/2) covariance function can lead to more reliable 
results than all others when applying GP to reconstruc¬ 
tions using D measurements. 

Following Refs, [l^, in which the detailed techni¬ 
cal description of GP can be found, we reconstruct the 
interaction between dark en erg y and dark matter using 
the SNIa Union 2.1 data set [^. Before that we will first 
show the reliability of the GP method. 


B. Mock data 

To demonstrate the ability of the GP method to distin¬ 
guish different models and recover the correct behaviors 
of the models, we create mock data sets of future SNIa 
according to the Dark Energy Survey (DES) for two 
fiducial models: the standard AGDM model and a toy de¬ 
caying vacuum model: pde = 3aH with w = —1. Here 
we set Dmo = 0.3 for both. 

For the AGDM model, it is straightforward to calcu¬ 
late the Hubble parameter H/Hq and then obtain the 
simulated data oi D{z) easily. Because there is no inter¬ 
action, we just want to check whether the reconstruction 
with this simulated data can recover q = 0. 

As for the decaying vacuum model, on the one hand we 
should calculate the Hubble parameter and on the other 
hand, we must get q{z) for this fiducial model. The Hub¬ 
ble parameter is simply iJ/i7o = ^ — + + 

The fiducial interaction is qHo^ = —3(1 — Dmo)(l + 
z)H'H. The rest of the task is to reconstruct q{z) from 
the Hubble parameter, and then test how well the re¬ 
constructed q{z) agrees with the fiducial one. If the GP 
method can recover both of the fiducial models and has 
the ability to distinguish them, we can demonstrate that 
GP is a valid method in the reconstruction for our pro¬ 
pose. In the following we will create the mock data sets. 

The DES is expected to obtain high quality light curves 
for about 4000 SNe la from z = 0.05 to z = 1.2 in the 
next five years. From Table 14 in we can calcu¬ 
late the errors of D: ajj and the corresponding num¬ 
bers of SNe la for each redshift bin. At every redshift 


point z, D{z)sim is sampled from the normal distribution 
D{z)sim ^ N{D(z)fid,crD) where D{z)sim is the simu¬ 
lated data of D{z) and D{z)fid is the theoretical value 
from the fiducial model. Thus we create the mock data 
sets for these two fiducial model. For each of them we re¬ 
construct D{z) and its derivatives and then apply to the 
reconstruction of the interaction term q{z). Note that 
each of the reconstructed D{z) and its derivatives at ev¬ 
ery redshift point is a distribution with a mean value and 
the error regions. However, D(z) and its derivatives are 
not independent but correlated by a covariance between 
them. Using the reconstructed D{z) and its derivatives 
also the covariance matrix between them, we can apply 
Monte Garlo sampling to determine the q{z) in Eq. (I7|) 
for a given w at each point z which we want to recon¬ 
struct. The detailed description of the covariance matrix 
can also be found in (2^. 

We can see from Fig. [T]that the D{z) and its deriva¬ 
tives are reconstructed very well from the mock data sets 
assuming the AGDM model. The dashed blue line is the 
mean of the reconstruction and the shaded blue regions 
are the 68% and 95% G.L. of the reconstruction. The er¬ 
rors for higher derivative of D{z) are a litter larger than 
the lower ones. As expected, at higher redshifts the er¬ 
rors become large due to the poor quality data in that 
region. 

Since our reconstruction involves D"'[z) whose errors 
are even lager in the higher redshift, and will definitely 
lead to somehow uncontrollable large errors in high red¬ 
shift regions when apply it to the Monte Garlo sampling 
in Eq. ©• For a better show, we introduce a prefac¬ 
tor (1 -I- z)“” to g(z), that is, q{z) = g(z)(l -I- z)“". The 
choice of n is somehow arbitrary, we take n = 6 here. The 
motivation to introduce the prefactor is comprehensible: 
we just focus on the interaction in the low and medium 
redshift range because the quality of observed data in 
the higher redshift range is so poor that it gives a weak 
constraint on the reconstructed interaction. Such a pre¬ 
factor is just considered as a scale transformation with 
respect to redshift, which does not influence the recon¬ 
struction of q{z) in the low and medium redshift ranges 
significantly and provides a better show in the higher red¬ 
shift range. Moreover, our aim is to examine the evidence 
of the interaction, namely, we are testing the equality of 
the quantity with zero. As a consequence, we are free 
to do this without loss of generality like what Ref. [30] 
has done for the null test of the AGDM model. In the 
rest of this paper, we use q as our interaction term to 
test the interaction and we mainly focus on the low and 
medium redshift ranges because of the poor quality of 
the observed data in high redshift ranges which will not 
give a good constraint for the interaction term. 

Figured shows that the reconstructed interaction q{z) 
is consistent with the fiducial AGDM model nicely, falling 
in the Icr limit. The reconstruction in the redshift range 
from 0 to 0.6 is better than that in the high redshift where 
the error is large. Further, we see from Figs O and |3| 
that the reconstructions of the decaying vacuum model 
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FIG. 1: Gaussian processes reconstruction of D{z), D'{z) {top), and D”{z), D'"{z) (bottom) obtained from a mock data set of 
future DES and assuming the ACDM model with flmo = 0.3 (red line). The dashed blue line is the mean of the reconstruction 
and the shaded blue regions are the 68% and 95% C.L. of the reconstruction, respectively. 
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FIG. 2: Reconstruction of q{z) (dashed line) from the mock data set of future DES and assuming the ACDM model with 
flmo = 0.3. The shaded blue regions are the 68% and 95% C.L. of the reconstruction. 



also recover the fiducial model very well, and obviously 
deviate from the ACDM model. This shows that the 
Gaussian processes can capture both the two models very 
well and correctly distinguish between them. 


C. Reconstruction using Union 2.1 data 

We now appD the GP method to the real data, Union 
2.1 data sets [^, which contains 580 SNela data. Let us 


consider 


m — M + 5 log 



-25 = 51og[(l + z)D], 


( 10 ) 


with Hq = 70 km/(s Mpc), following [s^. Actually our 
results are not sensitive to the values of Hq. The values 
of D depend on both Hq and the absolute magnitude M. 
We can fix the Hq and only consider the uncertainties in 
M as adopted in Ref. [^. We transform the distance 
modulus m — M given in the data set to D and set the 
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FIG. 3: Gaussian processes reconstruction of D{z), D'{z) {top), and D”{z), D'"{z) (bottom) obtained from a mock data set of 
future DES and assuming a toy decaying vacuum model: pde = 3aH with w = —1 and fl„i 0 = 0.3 (green line). The dashed 
blue line is the mean of the reconstructions and the shaded blue regions are the 68% and 95% C.L. of the reconstruction, 
respectively. The ACDM model is also shown (red dotted). 
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FIG. 4: Reconstruction of q{z) from a mock data set of future DES and assuming a decaying-vacuum model: pde = 3aH with 
w = —1 and fimo = 0.3 (green line). The shaded blue regions are the 68% and 95% C.L. of the reconstruction. The ACDM 
model is also shown (red dotted). 


theoretical initial conditions D{z = 0) = 0 and D'{z = 
0) = 1. First, let us consider the decaying vacuum energy 
case with re = — I [^. Figure [5] shows the reconstructed 
D{z) and its derivatives from Union 2.1 data. The error 
is a little larger than that for the reconstruction of the 
DES survey due to the smaller number of SNe la and 
larger measurement errors. The reconstruction of the 
interaction q{z) is shown in Fig. [6l We see that both the 


distance D{z) and interaction q{z) are consistent with the 
ACDM model, which implies that there is no evidence for 
the existence of the interaction. 

We now change the equation of state w to see how 
much the differences of w will influence the output of 
our reconstruction. We consider the luCDM model with 
w = —1.006 ± 0.045 from the Planck 2015 [s^, the 
CPL parametrization w{a) = wq + rca(l — a) with wq = 
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FIG. 5: Gaussian precesses reconstruction of D{z), D'{z) (top), and D"{z), D"'{z) (bottom) obtained from Union 2.1 data sets. 
The shaded blue regions are the 68% and 95% G.L. of the reconstrnction. The ACDM model (red line) is also shown. 
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FIG. 6: Reconstrnction of q{z) from Union 2.1 data sets. The shaded blue regions are the 68% and 95% G.L. of the reconstruc¬ 
tion. The red line corresponds to the AGDM model. 


— 1.046^0.170 Wa = 0.14^0,76 from HST Cluster Su¬ 
pernova Survey 2011 and other constant equation 
of state: w = —0.7, —0.8, —0.9, —1.1, —1.2 and —1.3, 
respectively, for a comparison. All of the results are pre¬ 
sented in Fig. [3 Here we emphasize that we have con¬ 
sidered the effects of the errors of w on the reconstructed 
q{z) in the cases of (d) and (e) in Fig. [T] Since the con¬ 
stant w for the icCDM model from Planck 2015 is very 
close to —1, so the reconstruction is almost the same as 
the w = —1 case. While for the CPL case, the error of 
the reconstruction is a little larger than those in other 


cases because of the poor constraints on wq and Wa ■ We 
see that if w lies between —0.9 and —1.1, q{z) = 0 is 
captured within the 95% confidence region of the recon¬ 
struction, while the interaction is shown up if the equa¬ 
tion of state for dark energy deviates signihcantly from 
— 1. This shows the fact that there is a degeneracy be¬ 
tween the interaction and the equation of state of dark 
energy, as indicated in Eq. 0. 
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(e) w(a) = wo + Wa(l — a) 




FIG. 7: Reconstructions of g(z) for the wCDM model and CPL parametrization. (a) w = —0.7; (b) w = —0.8; (c) w = —0.9; 
(d) w = — 1.006±0.045 from Planck 2015; (e) iu(a) = wo + Wa(l — a) with wo = —1.0461q j 7 q and Wa = 0.14^Q ®g; (f) w = —1.1; 
(g) w = —1.2; (h) w = —1.3; The shaded blue regions are the 68% and 95% C.L. for the reconstruction. 







































IV. DISCUSSIONS AND CONCLUSIONS 

We have presented an approach to reconstruct the in¬ 
teraction between dark energy and dark matter by using 
Gaussian processes. To check the reliability of the GP 
method, we create mock data for two fiducial models, 
one is the AGDM model which has no interaction be¬ 
tween dark matter and dark energy, the other is a toy 
decaying vacuum energy model with pde = SaH, which 
is of the interaction. It shows that the reconstruction 
method by using the Gaussian process works well and 
can capture the features of these two models. 

We then applied the method to the real data from 
Union 2.1 data sets and reconstructed the distance and 
its derivatives, and then the interaction. It was found 
that for the decaying vacuum energy model with w = 
— 1 [l^, there is no evidence for the existence of the in¬ 
teraction, namely, the ACDM model is consistent with 
the Union 2.1 data sets within Icr limits. 

To check the influence of the equation of state of dark 
energy on the method, we have also considered several 
constant values of w from —0.7 to —1.3. The results 
show that q = 0 falls in 95% G.L. of the reconstruction if 
w lies between —0.9 and —1.1. The wCDM model with 
w = —1.006 ± 0.045 from Planck 2015 falls in this range 
and g = 0 is within the Icr limits. This also holds for the 
GPL parametrization case. However, as we can see from 


Figure [3 if w deviates obviously from —1, the interac¬ 
tion exists beyond 2a C.L.. This reflects the degeneracy 
between the interaction and the equation of state of dark 
energy. 

Note that in our reconstruction method, only the ob¬ 
servational data on the expansion history of the universe 
can be used. For example, some measurements of Hubble 
parameter can also be combined into the above recon¬ 
struction. It is certainly of interest to develop a recon¬ 
struction method for the interaction between dark mat¬ 
ter and dark energy, in which some kinetic data of the 
universe, for example, the growth factor of large scale 
structure, can be used for this aim. Our paper shows 
that the existence of the interaction between dark mat¬ 
ter and dark energy found in [l^ is mainly due to the 
data of redshift space distortion. 
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